Peierls ground state and excitations in the electron-lattice correlated system (EDO-TTF)2X 
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We investigate the exotic Peierls state in the one-dimensional organic compound (ED0-TTF)2X, wherein 
the Peierls transition is accompanied by the bending of molecules and also by a fourfold periodic array of 
charge disproportionation along the one-dimensional chain. Such a Peierls state, wherein the interplay between 
the electron correlation and the electron-phonon interaction takes an important role, is examined based on an 
extended Peierls— Holstein— Hubbard model that includes the alternation of the elastic energies for both the 
lattice distortion and the molecular deformation. The model reproduces the experimentally observed pattern of 
the charge disproportionation and there exists a metastable state wherein the energy takes a local minimum with 
respect to the lattice distortion and/or molecular deformation. Furthermore, we investigate the excited states 
for both the Peierls ground state and the metastable state by considering the soliton formation of electrons. It 
is shown that the soliton excitation from the metastable state costs energy that is much smaller than that of the 
Peierls state, where the former is followed only by the charge degree of freedom and the latter is followed by that 
of spin and charge. Based on these results, we discuss the exotic photoinduced phase found in (EDO-TTF)2PF6. 

PACS numbers: 71.10.Fd, 71.10.Hf, 71.10.Pm, 71.30.-(-h 
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I. INTRODUCTION 

The quasi-one-dimensional molecular compound, (EDO- 
TTF)2PF6, (Refs. (where EDO-TTF denotes 

ethylenedioxytetrathiafulvalene) has been studied as one of 
the central topics showing a photoinduced phase transitioui^ 
The 1/4-filled (EDO-TTF)2PF6 system exhibits a Peierls insu- 
lating state, which is accompanied by a charge disproportion- 
ation and the bending of the EDO-TTF molecules.'* Several 
experiments indicate that this Peierls transition comes from 
the cooperation effect of the molecular deformation, charge 
ordering, and anion ordering. The charge disproportionation 
along the one-dimensional chain exhibits a fourfold periodic- 
ity given by a periodic array of [0, 1 , 1,0]; i.e., the pattern is an 
alternation of two charge-rich sites and two charge-poor sites. 
A remarkable feature of the insulating state is the large bend- 
ing of the EDO-TTF molecules, which does exist even for 
the neutral single molecule. The degree of bending depends 
on the valence of the EDO-TTF molecule. Each molecule 
in the high-temperature metallic phase has a valence of +0.5, 
while the valences in the insulating state are estimated as +0.9 
and +0.1 for hole-rich and hole-poor sites, respectively.^ In 
the insulating state, the degree of bending is enhanced at the 
electron-rich site and is suppressed at the electron-poor site; 
i.e., the bending (flattening) of the molecule is observed in 
the electron-rich (-poor) site. However, the mechanism be- 
yond the kinetic-energy gain is requked in order to explain 
the bending and/or flattening of the molecules. 

Recent experimental studies have focused on the phase 
transition induced by a weak laser pulse; i.e., the photoin- 
duced phase transition, and suggest that (EDO-TTF)2PF6 ex- 
hibits a gigantic photoresponse in the low-temperature phascji 
In the photoinduced phase, the lattice distortion and the 
molecular deformation are relaxed and the system shows a 
metallic behavior. In addition, it has been discussed that the 
photoinduced phase shows an exotic behavior, and is different 
from that in the high-temperature metal phase.- From these 
experiments, it has been argued that a new phase, which is not 



possible as a ground state and is related to a metastable state, 
is achieved by the photoexcitation. Therefore it is of particular 
interest to investigate the mechanism of a metallic behavior in 
such a metastable state, which may be related to the electrotiic 
correlation. 

A theoretical investigation on the Peierls state in the one- 
dimensional (ID) quarter-filled system has been performed 
based on the Peierls— Hubbard model including the effect of 
the onsite Coulomb repulsion U and the intersite one ]/jiL2ii2 
It has been clarified that the spatial variation of the conven- 
tional charge-density wave at the quarter-filled system is of 
the site-centered type; i.e., the charge density takes a max- 
imum at a site. On the other hand, in the Peierls state of 
(EDO-TTF)2PFg, it takes a maximum at the location between 
two neighboring sites; i.e., it is of the bond-centered type. It 
has been shown, within the mean-field theory, that the ex- 
perimentally observed pattern of charge disproportionation 
can be reproduced by taking into account the alternation of 
the elastic energies." It has also been shown that, due to 
the effect of the elastic-energy modulation, the Peierls state 
and the charge-ordered (CO) state with a twofold periodic- 
ity compete with each other, and there exists the metastable 
state wherein the energy takes a local minimum with respect 
to the lattice distortion." The trigger of such a metastable 
state comes from the competition between the electron cor- 
relation and the electron-phonon interaction. The former fa- 
vors the CO state with a twofold periodicity, while the latter 
supports the charge disproportionation with a fourfold period- 
icity. In the present paper, we study both the Peierls ground 
state and the metastable state, in terms of the phase represen- 
tation based on bosonization. Furthermore, we also exam- 
ine the two kinds of excitations: (i) the excitation followed 
by the lattice relaxation and (ii) the purely electronic excita- 
tions. In §2, our model consisting of an electron-lattice system 
is given and a representation based on bosonization is intro- 
duced. The Peierls distortion with the bond order is taken by 
considering the alternation of the elastic constant for the lat- 
tice distortion. In §3, the ground state is analyzed in terms 
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FIG. 1: (Color online) (a) Crystal structure of (EDO- 
TTF)2PF6 at RT and 260 K. (b) The HOMO energy 
level of the EDO-TTF molecules as a function of the 
distortion parameter 5 (see text). 5 = corresponds to 
the case at RT and 5 = 1 to that at 260 K. (c) The over- 
lap integral as a function of 5. At RT, all molecules are 
equivalent and S1=S3, while Sl7^S2. 



of the phase variables of spin and charge to show the condi- 
tion for the existence of the metastable state. Using an ex- 
tended Peierls—Holstein— Hubbard model treated in terms of 
the bosonization method, we demonstrate how the first-order 
transition occurs where the result is qualitative the same as 
the previous mean-field result. ' ' In §4, the excitations from 
both the Peierls state and the metastable state are examined 
by calculating the soliton formation energy. It is shown that 
the latter is much smaller than that of the former. Section V is 
devoted to the summary and discussions. 

II. MODEL FOR THE EDO-TTF COMPOUND 

As a system exhibiting the correlated Peierls state, we 
consider a tight-binding Hubbard model coupled with lattice 
distortion and molecular deformation, through the electron- 
phonon interaction. In this section, by focusing on the Peierls 
state in (EDO-TTF)2PF6, we first examine the characteristic 
features of the molecular deformation and the lattice distortion 
of the EDO-TTF molecules, based on a quantum-chemical 
calculation. Next, we derive the model Hamiltonian of the ID 
quarter-filled Peierls— Holstein— Hubbard model, and we in- 
troduce the "phases" of the lattice distortion and of the molec- 
ular deformation. 



A. Molecular deformation in (EDO-TTF)2PF6 

It has been confirmed that (EDO-TTF)2PF6 exhibits a first- 
order transition at around 280 K,i The overlap integrals at 
room temperature (RT) and at 260 K have been estimated by 
an extended Hiickel calculation'^ based on an x-ray structure 
analysis. At room temperature, there is a very weak dimer- 
ization and the overlap integral is almost uniform along the 
stacking direction. However, at 260 K, there is a strong vari- 
ation among the overlap integrals, given by 51, 52, and 53, 
along the stacking direction,'* 

In order to make the situation clearer, we estimate the en- 
ergy levels of the highest occupied molecular orbitals (HO- 
MOs) for the flat and bent molecules, and also estimate the 
overlap integrals, as a function of the lattice distortion and/or 
molecular deformation. Here we simply perform the inter- 
polation of the coordinates for the respective atoms in the 
EDO-TTF molecules between RT and 260 K by using the rela- 
tion (x, y, z) = {x, y, 2;)rt + (Aa;, Ay, Az)S, where 6 is the 



distortion parameter and (Aa;, Ay, Az) = (x, y, 2)260 k — 
(x, y, z)rt. For both cases of RT and 260 K, the y axis is 
chosen along the stacking direction [the vertical axis in Fig. 
[lla)] and the x axis is set to the projection of the molecular- 
long direction to perpendicular to y [the horizontal axis in Fig. 
[Ua)]. At 5 = (1), the atomic coordinates at RT (260 K) 
are reproduced. Using this relation and the extended Hiickel 
method,'^ we examine the S dependence of the HOMO energy 
levels [Fig.[nb)] and the overlap integral [Fig.lTJc)]. Here we 
neglect the flipping disorder of the terminal ethylene group^- 
and use only the coordinates of the atoms with a high occupa- 
tion rate. We found that the HOMO energy level for the bent 
EDO-TTF molecule is sufficiently lowered due to the molec- 
ular deformation, and that the overlap integrals do not change 
their sign as a function of S. We estimate the transfer inte- 
grals from the overlap integrals by using the empirical relation 
ti — SiX 10 eV.— At 260 K, the transfer integrals are esti- 
mated as ti = 442 meV, ^2 = 230 meV, and = 131 meV. 
A notable feature of this compound is that the electron-rich 
EDO-TTF molecule tends to bend and the hole-rich EDO- 
TTF molecule tends to be flat [Figs. 2(b) and 2(c)],'' where 
the valence of the former (latter) is 0.1 (+0.9). The energy 
level of the HOMO is ep = -9.222 eV for the flat EDO-TTF 
molecule, while eb — —9.442 eV for the bent molecule. This 
lattice distortion pattern shows a strong tetramerization ut in 
addition to the lattice dimerization Ud, which are estimated 
from 

h = t + gput + gpud, (2.1a) 
t2 ^ gpud, (2.1b) 
h = t - gput + gpUd, (2.1c) 

where gp is the electron-phonon coupling constant. The lat- 
tice distortion parameter ut characterizes a modulation of 
fourfold periodicity; i.e., the lattice tetramerization, and the 
parameter Ud characterizes a modulation of twofold period- 
icity; i.e., the lattice dimerization. It can be found that the 
strength of the lattice tetramerization is gput = 156 meV, 
while that of the lattice dimerization is gpUd — 29 meV. 

We note that the spatial variation of the lattice distortion 
coiTesponds to neither the conventional 2kp charge-density 
wave^'^ state nor the dimer-Mott-nspin-Peierls (SP) stateriiS 
in which the lattice tetramerization takes a maximum at ti 
(i.e., gpUt > 0), and the lattice dimerization has a maximum 
at t2 (i.e., gpUd < 0). In the dimer-Mott-nSP state, which 
occurs in the presence of the lattice dimerization, the system 
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FIG. 2: (Color online) (a) The HOMO energy dif- 
ference Ae and the magnitude of lattice tetrameriza- 
tion and dimerization as a function of the distortion pa- 
rameter S. (b) Schematic view of crystal structure of 
(EDO-TTF)2PF6 at 260 K and (c) the corresponding en- 
ergy diagram, where Ae — 156 meV, ti — 443 meV, 
t2 = 230 meV, and ta = 131 meV. 



becomes effectively half filling and each electron is localized 
at the location of each lattice dimer. When the 2k p Peierls in- 
stability is taken into account, such a paramagnetic insulator 
changes into the non-magnetic SP state. Then, one finds 
gput > and gpUd < in the dimer-Mott + SP state. How- 
ever, in the case of (EDO-TTF)2PF6, both the lattice tetmmer- 
ization and the lattice dimerization take the maximum at ti 
(i.e., gpUt > and gpUd > 0). This lattice distortion pat- 
tern in (EDO-TTF)2PF6 is qualitatively different from that 
in the dimer-Mott + SP state and has not been discovered in 
previous theoretical studies. In this sense, the Peierls ground 
state of (EDO-TTF)2PF6 is quite exotic and, is examined, in 
the present paper, by constructing a model that can reproduce 
such a ground state. 



B. Model Hamiltonian 

Based on the above consideration, we introduce a ID ex- 
tended Hubbard model coupled with the lattice through the 
electron-phonon interaction. Since there are two molecules 
within the unit cell, we introduce two sites referred to / = 1 
and Z = 2. In this study, we are based on the hole picture; i.e., 
there is one carrier per two sites. Our Hamiltonian is given by 



H — Hc + Hc-ph + -ffph, (2.2) 
where the purely electronic part is {j ^ 1, . . . , N) 

(4,2,sCj + l,l,s +H.C.) 

3 

(2.3a) 



the electron-phonon coupling term iJo-pti is 

3,s 

+ ^9P{UJ,2 - (c]2,sCj + i,i^s + H.c.) 

3-s 

+ Yl 9HVj,i Uj j - i j , (2.3b) 
and the phonon term iJph is 

3 

, sr^, 2 , 2 \ 

+ — Z^Kl+^J-,2) 

3 

+ ^i^(uj,i,Uj,2,Wj,i,t;j,2)- (2.3c) 

3 

The parameter Vj^i represents the degree of molecular defor- 
mation wherein the shape of the molecule depends on its sign; 
i.e., the molecule is bent for Vjj > 0, while the molecule is 
flat for Vjj < 0. In Eq. ( |2.3a| l, Cj^i.s denotes the annihilation 
operator of an electron with spin s(=t, j) at the Ith site within 
the jth unit cell, and n^,/ ^ — c^j j ^Cj ; where t is the trans- 
fer energy and /i is the chemical potential. The parameters 
U{> 0) and V{> 0) denote the magnitudes for the on-site 
and nearest neighbor interactions. In Eq. ( I2.3bb . gn is the 
electron-phonon coupling constant of the on-site (Holstein) 
type, and the Peierls-type electron-phonon coupling constant 
is given by gp. In Eq. ( I2.3cl ). the parameters Kh and Kp 
are the conventional elastic constant for the molecular de- 
formation and the lattice distortion, respectively. The quan- 
tity F{uj^i,Uj,2,Vj^i,Vj^2) represents the specific features in 
(EDO-TTF)2PF6, which wifl be defined in Sec. II C. 

Here, we note that the present Hamiltonian can reasonably 
account for the relation between the valence of the molecule 
and the molecular deformation Vjj. For simplicity, we con- 
sider the case of the single molecule, wherein the Hamiltonian 
is given by 

Hi = gnv - + ^^HV^- (2.4) 
By minimizing Hi with respect to v, we find that the lattice 
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distortion depends on the charge on the molecule, e.g., 

9H 



2Kh 

9H 

2Kh 



for n — 



for n = 1. 



(2.5) 



Thus, we obtain that the neutral molecule (n = 0) becomes 
bent (v < 0); on the other hand, the ionic molecule (n — 1) 
becomes flat (v > 0). 



C. Modifications of phonon term //, 



ph 



Here we examine the characteristic features of the crystal 
structure of (EDO-TTF)2PF6 [Fig.Ea)], for the model that re- 
produces the Peierls ground state. If F{uj^i,Uj,2, f j.2) = 
in Eq. (|2.3ct , the model reduces to the conventional 
Peierls— Holstein— Hubbard model at quarter filling .^i^i^i'° 
Noting two kinds of molecules in the unit cell at the high- 
temperature phase of the EDO-TTF compounds, we introduce 
two kinds of elastic constants for the lattice distortion and 
the molecule deformation, within the unit-cell and between 
neighboring cells. Such a difference would play important 
roles on the low-energy properties. We consider the term. 



■',2,Vj^i,Vj^2), given by 



F(Uj,l,Uj-,2,Wj,l,Wj,2) 

dKp 



5K 



(2.6) 



where the 6Kp term comes from the difference between the 
elastic constants for the lattice distortion, within the unit cell 
and between neighboring cells. The 5Kh term represents the 
difference in the energy gain of the two-molecular deforma- 
tion within the unit cell. The intrinsic property of the bending 
of the EDO-TTF molecules would play important roles for 
the lattice distortion, since a pair of two molecules faces each 
other in the unit cell.'* In fact, in the ordered phase, one pair 
of molecules shows strong bending and another pair of the 
molecules becomes more flat [Fig. [Ha)]. Such an effect due 
to the pairing of the EDO-TTF molecules can be taken into 
account in our model by considering the additional terms of 
Eq. (ESll. 

From the mean-field approach," it has been pointed out 
that the effect of 5Kp is important to determine the multi- 
stable states induced by the lattice distortion. The details of 
the effects of 5Kp and 5Kh are discussed in Sec. III. 



D. Lattice distortion and molecular deformation 

In general, the lattice distortion Uj_i in quarter-filled sys- 
tems can be decomposed into the component with a twofold 
periodicity (i.e., dimerization), u^, and that with a fourfold pe- 
riodicity (i.e., tetramerization), ut-'^''^'^^ Even for the freedom 



of the molecular deformation, Vj,i, it can be decomposed into 
the component with a twofold periodicity, v^, and that with 
a fourfold periodicity, vt- We introduce these quantities by 
using the relations: 



Ut 



+ (-l)^^cosC- 



Ud 

2 ' 



+ (-lF^COsC-y, 



(2.7a) 
(2.7b) 
(2.7c) 
(2.7d) 



The phases ^ and ( determine the spatial pattern of the 
tetramerizations ut and vt, respectively. 

For the later convenience, we renumber the site index j as 
i = 2j -\- {I — 1); accordingly, we rewrite Cjd.s Ci,s- By 
using Eq. (I2.7l i. the Hamiltonian ( I2.2l l is rewritten as 



H 



Y[t- gpUt cos (I* + C - ^) + (-l)'5P^id 



4 * 



SKp n Kp n 

N^u'^sin2^ + N^ul 



N 



Kh - SKh 2 ixT^^H 2 
: Vf ~ Al 



^ v^sin2^ + N^vl 
(2.8) 



Now, we specify the phases of ^ and (, which repro- 
duce the Peierls state observed in (EDO-TTF)2PFe. The 
crystal structure of (EDO-TTF)2PF6 for the metallic state 
at high temperature shows an alternation of the bending 
of the molecule, where all of the molecules are identical 
crystallographically.^ For the insulating state at low temper- 
ature, the neutral molecule exhibits a large bending and the 
ionic molecule becomes rather flat, where the adjacent two 
molecules at the neutral sites are oppositely located with the 
same degree of bending and becomes convex outside. The 
transfer integral between adjacent bent molecules becomes 
small while that between adjacent flat molecules becomes 
large. Such a pattern of the molecular deformation and lat- 
tice distortion can be reproduced by setting ^ — C, — tt/A. 
When ^ = C = tt/A, the transfer-energy modulation is given 
by 



(^3,^2,^1,^2,^3, ^2,^1, ^2), 



(2.9) 



for i = 0, 1, 2, 3, 4, . . ., where ti, t2, and are the same as 
those given in Eq. ( 12.1b . The schematic view of the modula- 
tion pattern is shown in Fig.|3] The site-energy modulation is 
given by 

7:9H(+vt, +vt, -Vt, -vt,-\-vt, +vt, -Vt, -vt, ■ ■■), (2.10) 
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for i = 0, 1, 2, 3,4,.. ., where the component of the molec- 
ular deformation with the twofold periodicity vanishes for 
(EDT-TTF)2PF6 (see also Fig.O. 

III. BOSONIZATION REPRESENTATION 

Here, we represent the Hamiltonian in terms of bosonic 
phase variables jii By setting a as the lattice constant, the elec- 
tron density operator can be represented as^^ 

Tij 11 d(j)„ 2 . 

— = \ s\n(2k F X ^- (pp) cos (pa 

a Aa n ax ira 

2c 

cosiAkpx + 2(l)p), (3.1) 

na 

where x = ia and c is a nonuniversal constant. The Fermi 
momentum is fc^ = 7r/(4a). The quantities cpp and (p^ are 
bosonic phase variables and Eq. ( 13.1b is justified by the mi- 
croscopic representation of the field operator 

Here we briefly recall the relation between the lattice dis- 
tortion and electron density modulation, based on the phase- 
variable representation. We tentatively assume c = 0. If the 
phases are locked at {{(fip), {<j>a-)) — i'^/^:, 0), the expectation 
value of the charge-density operator are given by 

1 2 

(n^) = sin(2fci^a; + 7r/4). (3.2) 

4 TT 



The pattern of charge modulation is given by (rij) = (o o 
00° • • ■ ) fori = 0, 1, 2, 3, 4, . . ., which is compatible with 
the lattice distortion and molecular deformation given in Eqs. 
i2.9i and ( 12.101 ). This means that the tetramerization pattern 
observed in (EDO-TTF)2PF6 is reproduced if the phases are 
locked at e = C = 7r/4and ((0p), (0,)) = (^4,0). 



Based on the phase-variable representation, the bosonized 
Hamiltoniani^ii^^ is obtained as H ~ J dxH, where 



Uf O O O O 0*D O— 





I Ki. 3: (Color online) The modulation pattern of the transfer inte- 
gral induced by the lattice distortion ut and that of the onsite energy 
level induced by the molecular deformation vt . In this view, we have 
considered only the modulations of tetramerization; i.e., of fourfold 
periodicity. 



Vp_ 
477 



■^{d^cl>pf + Kp{d^0pf 

Kp 



47r 

1 



1 



gi/4 

27r2a2 



cos 40p + ^ cos 2(j)a 



- —gput cos {(j)p - £,) cos (pa ^ — gnvt sin {(j)p + C) cos (f>a + -^gpu^ sin 2(j)p + —^——gnVd cos 2(p 

Kp 2 SKp 2 . r,c , 2 , Kh - SKh 2 ^I^H 2 ■ 

+ ^'t n — ^'t sm 2^ + — — H 7i — sm 2C, + -— ■ 

4a ' 8a ' 2a 8a 8a 8a 



27r2a2 
U 



Kh 2 



r 



(3.3) 



In Eq. ( 13.3b . 6^ (v = p and a) is the phase variable 
which is canonically conjugate to (f>i, and 6',y(x')] = 

i7rsgn(a:: — x'). The parameters Vp and Va are the velocity 
of the charge and spin excitations, and Kp and K^ are the 
Tomonaga—Luttinger parameters. It is noted that Kp < 1 for 
the repulsive interaction and Ka- — 1 for the paramagnetic 
state. From the perturbative calculation,'^ the coupling 51/4 
is given by 171/4 oc U^{U — 4y), where it is noticed that gi/4 
is positive for small V but becomes negative for large V. For 
small V, the phase (pp is locked al {(f>p) — 7r/4(mod 7r/2); 
while for large V the phase is locked at {(pp) — O(mod 7r/2). 



In terms of the bosonized Hamiltonian, we examine the ef- 
fect of additional terms including &Kp and SKh- First, we 
discuss the state in the case of 5Kp = SKh = 0. The lock- 
ing position of the phases ^ and ( is determined by (/)p; i.e.. 



and C = 7r/2 



Thus, the total Hamiltonian can 



be expressed as 



'H\sK=o ~ -r 

4:71 



— {d^ApY + Kp{d,6p) 
Kp 



— {d^(j)af + Kp{dJaf 



cos 4(/)p + 



in 

gi/4 

27r2a2 

— gpUt cos (pa 
na 

U 



da 



27r2a2 
1 

^/2'n 



cos 2(pa 



-QHVt cos (pa 



2 opurf sin 20P 



u 



TT-^ta 



2V2TT^ta 



gnvd cos 20P 



Kp 2 Kp 2 , Kh 2 , Kh 2 
4a 2a 8a 8a 

which can reproduce the previous results (e.g.. Fig. 4 in Ref. 
8) qualitatively. The undistorted state can be obtained for a 
small electron-phonon coupling of gp and gn- By increasing 
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gp only, we obtain a solution where Ud becomes finite and 
the phases (j)p and ^ are locked at (/)p = ^ = 7r/4(mod tt) (if 
Ud < 0) or at (/)p = ^ = 37r/4(mod tt) (if > 0). This 
state corresponds to the Akp bond-order wave (BOW) state in 
Ref. Bt. By further increasing gp, we obtain a state wherein 
the tetramerizations ut and vt become finite. For Ud < 0, the 
lattice distortion [Eq. (|2.7| i1 is given by 



2 ^ 2 ^ 



2 



(-1) 



j,2 2 ' ' 2 

Then the resultant hopping integrals are given by 



ti ^ t- gp\ud \ + gput, 

h = t + gp\ud\, 

h ^ t- 9p\ud\ - gpUt- 



(3.5a) 
(3.5b) 



(3.6a) 
(3.6b) 
(3.6c) 



where h ^ [t - gp{ui,i - ^1,2)], h = [t - 5p('"i,2 - "2,1)]' 
and = [t — gp{uj+i^i — Uj^i^2)]- We note that this state 
is nothing but the bond-charge-density wave (BCDW) state in 
Ref. 8 and the dimer-Mott-nSP state in Ref.|9t 

Next we examine the effect of 5Kp and 5Kh From Eq. 
(13. 3t . we can immediately find that, if 5Kp > 0, the phase 
^ tends to be locked at ^ = 7r/4(mod tt). In addition, the 
phase C, also tends to be locked at C, — 7r/4(mod tt) due to 
5Kh{> 0). From Eq. (I3.3l l. we can immediately find that the 
most plausible set of the locking pattern of the phases in the 
Peierls state is given by C, 0p, = (7r/4, 7r/4, 7r/4, 0). 
Thus, the corresponding charge modulation favors the bond- 
centered type. In the following analysis, we restrict ourselves 
to the case of ^ = C, — it /A. Another effect of 5Kp and 
5Kh is to decrease the elastic energy of the tetramerizations 
Ut and vu i.e., the tetramerization is favorable compared to 
the dimerization Ud and the alternation Vd- Thus we discard 
Ud and Vd in the following analysis. In this case, the total 
Hamiltonian can be simplified as 



n 



Att 



Att 

1 

2^a 



(3.7) 



where the potential term is given by 



V{ut,4>p,4>a) = .91/4 COS 40p + COS 2(/)cr 



Uf COS 



COS ( 



(3.8) 



Here we have rescaled {4:TTa)[gpUt + g/fWf /(2V2)] 
Ut and ifph is evaluated from the variational method as 
Kph = [{Kp + KH)/{KpKH)]/ma) where Kp = (Kp - 
5Kp/2)/gl and Kh = 4:{Kh - 25KH)/g%. Hereafter we 
perform the classical treatment for the qualitative understand- 
ing of the metastable state at zero temperature. 
IV. PEIERLS STATE VERSUS CHARGE-ORDERED STATE 



In this section, the ground state and the metastable state are 
examined by calculating the minimum energy as a function of 
Ut- Since the phase variables (pp, (f>a-, and the lattice distortion 
Ut are spatially uniform, we can focus only on the potential 
term V{ut, 4>p,(t>a) by discarding the first and second terms 
of Eq. ( 13.71 ). The ut dependence of the ground-state energy 
is estimated by using the stationary conditions for the phase 
variables, which are written as(^v — p and a) 



dV{utApA<y) 



0. 



(4.1) 



First we consider the case of V = Q. In this case, the 
coupling constant of the commensurability energy 31/4 be- 
comes positive, and then the charge phase 0p is locked at 
(/)pP' — tt/A. As for the spin phase (pa, the optimized locking 
position is 0°p' = cos~^ {ut/uD for Ut < , and (/)°p* — 
for Ut > u^, where uf = Ag„. Then, the minimum energy 
E^{ut) = yo(ut,CP\CP')/(2^'a') is given by 



Eo{ut) 




|gi/4| +9a-'' 



ph 2 

—u, 



for < ut < Ut 



for Uf < Ut- 



(4.2) 



From the coefficient of Uj of Eq. (I4.2l i. it is found that the state ut = becomes unstable and the Peierls state is obtained for 
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FIG. 4: The energy Eo{ut) as a function of ut, where itf — 'ig^ > 
0, and we have set Qi/i/ga = +1/4. The elastic constant is fixed as 
ffo-A'ph = 1/3, 1/4, 1/5, and 1/8, from top to bottom. The critical 
point is QaKpi, = 1/4. 



• 


Peierls 


o 


CO 



O 6 c> - - O 



o -o --0-- o- — o 



(t>a in the limits of the Peierls and CO states. The arrow indicates the 
schematic variation of the locking potential with increasing Ut . 



ga < l/{4:Kpi,). On the other hand, for > l/{4:Kpi,), the 
undistorted state (ut — 0) is obtained, leading to the spin- 
density-wave (SDW) state. The typical ut dependences of 
i?o(ut) are shown in Fig.|4] 

For the case of large V, the coupling constant of the com- 
mensurability energy gi/4 changes its sign (51/4 < 0), and 



favors the locking position (f>p = O(mod tt), which is noth- 
ing but the CO state with a twofold periodicity; i.e., {uj) — 
( O ° O ° ■ ■ ■ )■ III this case, the region for the lattice dis- 
tortion in divided by two kinds of characteristic values and 



which are given by 



\/QM9i/4\9<t, 



^9o 



For small ut{< u^,Uf), the phases (pp and 
at intermediate values, while the phase 0^ takes {(pp, 
(7r/4,0), (-37r/4,±7r), .... for large ut{> u^t^uf). For 
the intermediate value of ut, there are two possibilities: (i) 



(4.3) 

are locked 

) = 



< uf and (ii) u: 



on the former case, 



t < "t- 



In the present paper, we focus 
< (which is satisfied for gi /4 < 
and |gi/4| < (7cr/4), since within the mean-field calculation,— 
only the former case is realized and there would be some sub- 
tleties in the case of |gi/4| > 5cr/4. 

For case (i), the locking positions for the phase variables <j>p 
and (pa- are analytically determined as 



■ sm 



,P2 



TT 

4 



for < Ut < 
for u1 < Ut, 



(4.4) 



and 



-opt 




for Q < Ut < Ut 
for u1 < Ut < Ut 

\ "t / 

for Ut < Ut- 

(4.5) 

The schematic variations of (/)°p* and 0°p* as a function of Ut 
are shown in Fig.|5] 

Now, we examine the energy as a function of ut- By insert- 
ing Eqs. ( 14.41 ) and ( |4.5l l into Eq. ( 13.81 ), the energy with fixed 

Ut is given by Eo{ut) = Vo{ut, cP°p'\(P°P')/i2n^a^), where 



Eo{ut) = < 




\9l/4\9a 



for < Ut < Ut 



for Ut < Ut < Ut 
for Ut < Ut- 



(4.6) 



Typical ut dependencies of Eo{ut) are shown in Fig. |6l For ut > u^, the energy Ed(ut) takes a same form given in Fig.|4] 
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FIG. 6: The energy Eoiut) as a function of Ut, where Ut = 4(7^, 
and we have set gi/n/ga = —1/4, The elastic constant is fixed as 
QcrKph = 1/3, 1/4, 1/5, and 1/8, from top to bottom. The transition 
point is QirKph — 1/5. 
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It is noted that, in the present case, the energy around ut = 
can become a local minimum; i.e., the undistorted state is re- 
alized as a metastable state. The actual ground state is deter- 
mined in order to minimize the energy i?o(ut) with respect to 
Ut- From Eq. ( 14.61 ), we find that the Peierls state with finite ut 
is obtained in the case of 



9t + I51/4I < 



4Kph 

The amplitude of the lattice distortion is given by 

1 



ph 



(4.7) 



(4.8) 



In addition, the condition for the Peierls state with the 
metastable CO state (ut = Oj is given by 



I51/4I < 



1 



ph 



< 2.9^- 



(4.9) 



By further increasing Kp\^{ > 1/ [4(5cr + |gi/4 1 )] }, we find that 
the energy at = becomes lower than that of finite ut \ i.e., 
the first-order phase transition from the Peierls state into the 
CO state occurs when l/{4,Kph) = + \gi/i\. 

We note that the potential barriers from the Peierls state, 
AEp = Eo{u^) - £'o(^t°), and that from the CO state, 
A£^co = -£'o(Wf ) — i?o(0), are respectively given by 



AEp 



i9| 



AEco 



where ui 



2"|5i/4l5', 
1 



2I51/4I - 25^ 

i^ph 1 



1 



2i^ph 

2' 



165^ 



27r2a2 



ph 



(4.10a) 



, (4.10b) 



. 2'|ffi/4r/'5. V^ph/2 - 1/(16.9.) is the point 
at which the energy takes a local maximum. 

Based on the calculation of the ground state energy, we ob- 
tain the phase diagram on the plane of ^4 and Kpi^ shown 
in Fig.|7] For 51/4 > in which the CO state is absent, the 



FIG. 7: (color online). Ground-state phase diagram on the plane of 
51/4/ffCT andgcifph. In the shaded region, there appears a metastable 
state. 



pure SDW state is obtained for large Kp^ and the Peierls state 
is obtained for small Kph- When 171/4 < 0, the Peierls state 
is obtained for 4/\ph < l/{ga + 131/41) (region I) and the CO 
state coexisting with the SDW state is obtained for 4A'ph > 
l/(.9cr + I51/4I) (region II). The metastable state at itj = is 
obtained for iiTph > 1/ (S^cr) in region I, while the metastable 
state at U( 7^ is obtained for Xph < 1/(4(7.) in region 
II. By taking into account the relations 171/4 cx U'^{U — W) 
and ga oc Ur^ we can obtain the ground-state phase diagram 
on the plane of V and the elastic constant, which reproduces 
that of the mean-field theory. ' ' In section V, we examine the 
excited state in region I where the metastable state exists at 
Ut — 0, as shown in Fig.|6] 



V. EXCITATIONS BY SOLITON FORMATIONS 

In this section, we examine the excitations due to the soliton 
formations. First, we derive the equations for spatially depen- 
dent phases, which describe the variation from those of the 
Peierls ground state and the metastable state around ut — 0. 
As shown in section TV, the phase variables are uniform in 
space in the ground state. For simplicity, we do not consider 
the effect of the quantum fluctuations induced by 6p and 6'.. 
Since such effects give rise to a reduction in the amplitude 
of the potential V{ut,(f>p,<j)a-), it rnay be qualitatively under- 
stood by introducing the effective coupling constants given 



by 5?f4 = 5i/4exph8((/)2)], gf ^ 5, exp[-2(02)], and 



,,cff 



— Ucff exp[— i(((/)p) + (0^))], which can be estimated 
using the renormalization group treatment. '^ 

We examine the classical Hamiltonian given by Hc\ = 
J dxHci, where 



4ttK„ 



d(j)p{x) 
dx 



d4>a{x) 

dx 



-V{ut{x),(j>p{x),(l)a{x)). 



(5.1) 
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Minimizing Eq. ( 15. Il l with respect to ut{x), 4>p{x), and 
(f)„{x),^'' we obtain the following equations: 

ut{x) = -^cofi{^(t)p{x) -^cos4)„{x), (5.2a) 

° = ~t + 4.91/4 sin 40,(x) 

- ut{x) sin (^(/'^(a;) - cos (/)cr(a;), (5.2b) 
= — — — + 2ga- sm20<^(x) 

- ut{x) sin (^4>p{x) - sin0CT(a;), (5.2c) 

which are self-consistently determined. Note that substituting 
Eq. ( I5.2al i into Eqs. (I5.2bl i and (I5.2cl i. we obtain the equations 
for the phase variables (j)p{x) and (paix). 

In the conventional picture of the soliton excitations from 
the Peierls ground state, ''^ the charge and/or spin solitons are 
accompanied with the soliton of lattice distortion which con- 
nects the two minima of the potential. Due to this effect, the 
midgap state appears in the presence of the lattice distortion. 
In the present case, however, we have another excitation due 
to the lattice relaxation. As seen in Sec. IV, such an undis- 
torted state (ut — 0) can be locally stabilized due to the com- 
petition between the Peierls state and the CO state. In this 
section, first we verify that such a metastable state can be sta- 
bilized within the reasonable choice of parameters. Next, we 
also consider the purely electronic excitations by considering 
the antiadiabatic limit where the lattice distortion is assumed 
to be uniform, and we estimate the magnitude of the excita- 
tion gap. These states are calculated to comprehend the state 
relevant to reflectivity measurements in the Peierls phase and 
in the photoinduced phase in (EDO-TTF)2PFg. 

A. Lattice relaxations: Tlie appearance of the cliarge-ordered 
domain in the Peierls state 

We consider the case wherein both the Peierls state and the 
CO state become locally stable. Even for the case in which 
the true ground state is given by the Peierls state, it is actually 
expected that the Peierls state coexists with the metastable CO 
state without lattice distortion, in the photoexcited phase of 
(EDO-TTF)2PF6. 

Figure|8]shows the coexisting state between the Peierls state 
and the SDWh-CO state, created by the soliton-antisoliton for- 
mation connecting the different two states, which is obtained 
by solving Eq. iS.H . The parameters are set as gi/i/ga = 
—0.25 and gcr^^ph = 0.2, which correspond to the values 
on the boundary (see Fig. |7]i. The parameters for the ki- 
netic part are chosen as Kp = 0.1, K„ = 1, and Vp = 
Va — V^t. There exists a domain of the metastable CO 
state with {ut, (f>p, 4>a) = (0, 0, tt/2) in the Peierls state with 
{ut, 4>p, (f>a) = 7 ''■/4, 0). Note that the boundary between 
the Peierls and the CO states is created by the formation of 
solitons for charge, spin, and lattice distortion. Since such a 
coexisting state is expected not only for the Peierls phase at 




FIG. 8: (Color online) (a) The lattice distortion for the coexisting 
CO and Peierls states (the schematic view is also shown) and (b) the 
corresponding soliton-antisoliton formation in the phase variables (fip 
and We have set gi/4,/ga = —1/4, Kph — l/{5ga), Kp = 0.1, 
A'ct = 1, and Vp = Va- — %/2t with ~ 2t. 



finite temperature but also for the state after the photoinduced 
phase transition in (EDO-TTF)2PF6,' it can be expected that 
such a domain stays as long as the thermal fluctuations cannot 
go over the potential baiTier given by Eq. ( I4.10al ). 



B. Purely electronic excitations: the antiadiabatic limit 

Now, we examine the excited states induced only by elec- 
tronic excitations; i.e., spin and charge soliton-antisoliton ex- 
citations, while we assume that the lattice distortion remains 
uniform in space. This situation could be related to the elec- 
tronic excitations due to the probe light in reflectivity mea- 
surements. Such a soliton excitation is also expected to play 
crucial roles in the photoinduced phase transition in (EDO- 
TTF)2PF6 due to the ultrafast photoresponse within the order 
of pico seconds. In general,— the lattice should locally relax 
under friction to the minimum of the adiabatic potential, when 
photons are absorbed by electrons at a site. 

To this end, we calculate Eqs. ( |5.2b| ) and ( |5.2c| ) with the 
fixed uniform ut, for the Peierls state (ut ~ u^) and the CO 
state (ut — 0). In the Peierls state corresponding to finite ut 
(= the soliton-antisoliton excitation is characterized as 

(0P>.)=(J,O)^(-^^,.)^(5,O). (5.3) 

The explicit profile is shown in Fig. |9|a). Thus, the soliton- 
antisoliton excitation in the Peierls state carries the topologi- 
cal charge Q — ±1 and the spin Sz ~ ±1/2,^*^ which is noth- 
ing but a single-electron excitation. The simple picture is the 
following: A particle at the position of the soliton moves away 
from the place of the antisoliton by breaking the spin-singlet 
state. Then, one finds an emergence of a local spin with a hole 
around the soliton kink and that of a local spin with an extra 
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(a) Peierls state 




(b) CO state 



50 
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100 



FIG. 9: (Color online) Soliton-antisoliton excitation of the charge 
and spin phase variables, with fixed several values of (a) ut = m? 
(corresponding to the Peierls state) and (b) iit = (corresponding to 
the CO state). The soliton and/or antisoliton of <^p and <^o- connects 
two neighboring minima of V{ut, (l>p, (l>cr) (see Fig. [Sjl. The low- 
energy excitation in the Peierls state is a single-electron excitation, 
while that in the CO state is a domain excitation. The schematic 
views of the excitations for the respective cases are also shown. We 
have set gi/i/ga — —1/4 and A'ph = 3/(16(?o-) with ga — 2t, 
Kp = 0.1, and = 1. 



electron around the antisoliton kink. On the other hand, the 
soliton in the CO state is given by [Fig.|9jb)], 



TT vr 
2' 2 



(5.4) 



which carries the topological charge Q — ±1/2 and Sz — 0. 
A noticeable feature is seen from the fact that only the charge 
excitation occurs. This is a typical domain excitation in the 
CO state. The charge increases at the kink of the soliton and 
the extra hole appears at the kink of the antisoliton. 

As shown in Fig.|5] the locking location of the phase vari- 
ables (f)p and (f)„ changes depending on the strength of the lat- 
tice distortion ut- Even though the state with the intermediate 
values of < Ut < is not stable, we can consider the 
purely electronic soliton excitations in such a virtual state and 
discuss the "uf dependence" of the soliton excitation energy. 
The Ut dependence of the ground-state energy and the soliton- 
antisohton energy is shown in Fig.[lO] Here, we assumed that 
the soliton and antisoliton do not form a bound state (breather) 
and estimated the soliton-antisoliton formation energy by dou- 
bling the soliton energy in Fig. |9] From Fig. [TOl it is found 
that the soliton formation energy in the CO state costs less en- 
ergy compared with that in the Peierls state. This is due to 
the fact that only the charge degree of freedom participates 
in for the CO state. From this excitation energy analysis, it 




FIG. 10: (Color online) Ground-state {Eq) and soliton-antisoliton 
formation {Es) energies as a function of ut, where we have set 

gi/ilg<y = -1/4 and A'ph = 3/(16g<j) with g^ = 2t, Kp = 0.1, 
and Ka = 1. 



can be found that both the Peierls state {ut = u^) and the CO 
state {ut = 0) have finite energy gaps to the excited states; 
i.e., both states are insulators. However, in the CO state, the 
system could have a semiconducting feature since the gap can 
become very small. The magnitude of this energy gap depends 
on the choice of parameters, but it is discussed by comparing 
the mean-field results in Sec. I VII 

There are many choices of creating solitons in Fig.|9j how- 
ever, the present estimation has been done for the soliton, 
which has the minimum energy for large ut limit and is adia- 
batically connected to the small Ut case. From this excitation 
energy analysis, it can be found that the one-particle excita- 
tions in the Peierls state could be related to the domain excita- 
tions in the CO state. Such a connectivity would be relevant to 
the photoinduced phase transition mechanism, but this needs 
further investigation. 



VI. SUMMARY AND DISCUSSIONS 

We have examined the exotic Peierls state observed in 
(EDO-TTF)2PFg and analyzed the excitation from the Peierls 
ground state and also the excitation from the metastable CO 
state. The metastable CO state around the CO state has been 
found by considering the modulation of the elastic constants 
5Kp and 5Kh and by considering the moderate strength of 
the nearest neighbor repulsive interaction. The metastable 
state, which coexists with the spin-density wave, exhibits soli- 
ton excitations, having an energy much smaller than that of 
the Peierls state. The energy of the ground state and the ex- 
cited state, which are obtained as a function of ut, exhibits a 
common feature with the adiabatic potential in the electron- 
phonon coupled system in one dimensioni^ Considering the 
correspondence between such a situation to the case of EDO- 
TTF, the present metastable state could be relevant to the state 
of the photoinduced phase. 

Here, we discuss the experimental findings on the reflectiv- 
ity in the photoinduced phase of (EDO-TTF)2PFg based on 
the present analysis. From the photoconductivity measure- 
ment, the photoinduced phase exhibits a metallic behavior, al- 
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though it is not the same as that of high-temperature metallic 
state, as shown by the reflectivity data.^ In the present analy- 
sis, we obtained the metastable state without lattice distortion, 
which is not the normal state. Actually, in such a state, the 
translational invariance is broken due to the presence of the 
CO state, which gives a small soliton excitation gap. Here, we 
estimate the magnitude of the soliton gaps. From the mean- 
field approach,^^ the commensurability potential is estimated 
as Co cos40p, where the amplitude is given by Cq « — 0.03t 
for {U/t, V/t) = (4, 2). On the other hand, the amphtude of 
the cos2(/)ct potential^-' is about O.lSt ~ 0.3i, which are es- 
timated from {V/t,V2/t) K (1.4, 1.05) - (2, 0). Thus, we 
fmd gi/^l ga w —0.1 ^ —0.2 and our choice of the parameter 
9i/il 9<y = ^1/4 in Fig.[TO]is not far from the mean-field re- 
sults. Actually, we obtain ga/i^ir'^a^) — — 45i/4/(27r^a^) = 
Q.l2t from the mean-field result, Ca{V = 2t) « -0.031 
Here, we note the soliton gap at ut — 0, which is of the 
order of 5E{Q) = [£'soiiton(0) - i^o(O)] w O.Oli, and then 
5E{Q) Ri 3 meV. On the other hand, the soliton gap from the 
Peierls ground state is given by SE{u°^^) = [£'soiiton(u°'^*) — 
Eq{u°^^)] « 0.045t ~ 14 meV, which is much larger than 
that of the metastable state. 

Finally, we comment on the quantum fluctuation on the lat- 
tice distortion, which is not treated in the present calculation. 
Such an effect appears when the term J2j Pj /2M is taken 
into account for Eq. ( |2.3ct with [uj, Pji] = Sjji. The present 
calculation corresponds to the limit of a large mass M. If 
we naively consider a renormalization by a Debye— Waller 



factor, the effective coupling constant is reduced due 
to the quantum fluctuation. Thus, this fluctuation makes the 
normalized constant it^ph large and the region for CO in Fig. 
|7]is extended. For the limit of the light mass, which corre- 
sponds to the antiadiabatic case, the averaged lattice distortion 
vanishes but the crossover from the spin-density-wave state to 
the spin singlet is expected with an increase in the electron- 
phonon coupling constant .gph(> \/KphU) -^^'^^ Thus, the 
metastable state, which coexists with the CO and spin-density- 
wave states due to the electronic correlations, is expected to 
give rise to the conduction by forming the solitary charge ex- 
citation if the quantum fluctuations of both the electron and 
the phonon remain moderately small. 
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